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Using the finite-element discrete variable representation of the nonequilibrium Green's function (NEGF) 
we extend previous work [K. Balzer et al., Phys. Rev. A 81, 022510 (2010)] to nonequilibrium situations 
and compute — from the two-time Schwinger-Keldysh-Kadanoff-Baym equations — the response of the helium 
atom and the heteronuclear molecule lithium hydride to laser fields in the uv and xuv regime. In particular, by 
comparing the one-electron density and the dipole moment to time-dependent Hartree-Fock results on one hand 
and the full solution of the time-dependent Schrodinger equation on the other hand, we demonstrate that the 
time-dependent second Born approximation carries valuable information about electron-electron correlation 
effects. Also, we outline an efficient distributed memory concept which enables a parallel and well scalable 
algorithm for computing the NEGF in the two-time domain. 
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I. INTRODUCTION 

Owing to binary interactions, electrons inside and between 
atoms and molecules form a highly sophisticated quantum 
many-particle system. Already two-electron atoms show rich 
structure (T), e.g. revealed by the presence of autoionization 
resonances above the lowest ionization threshold. Therefore, 
when out of equilibrium, and, especially, when exposed to 
intense laser fields, atoms and molecules comprise very com- 
plex electron dynamics which occur at sub-femtosecond time 
scales and include nonlinear excitations, electron scattering 
and ionization fl2). Novel techniques, such as attosecond streak- 
ing (using coherent, ultrashort laser pulses) [3 1, allow one to 
capture and trace this dynamics experimentally and give, e.g., 
time- and angle-resolved indication of single, double or multi 
ionization processes. One of the most prominent examples, in 
this context, are attosecond (two-color) pump-probe scenarios, 
e.g. Eli including ionization that follows an electron shake 
up J5). Thereby, laser sources allow for the preparation of 
specially selected atomic or molecular initial states. 

In all aforementioned processes, electron-electron (e-e) cor- 
relations play a substantial role. Thus, consistent methods 
beyond the single-active electron picture are required to de- 
scribe electronic structure and dynamics. However, theory 
largely depends on good approximate solutions to the many- 
electron problem and, particularly, time-resolved investigations 
including correlations in sub-femtosecond regimes are still in 
their infancy. 

In this paper, we, in an ab-initio fashion, apply generalized 
kinetic equations, the Schwinger-Keldysh-Kadanoff-Baym 
equations (SKKBE) [6-8 1, to compute the electron motion 
in model atoms and molecules exposed to laser fields. In 
contrast to time-dependent density functional theory, NEGFs, 
thereby, offer a systematic approach to e-e correlations, where 
the residual interactions beyond the time-dependent Hartree- 
Fock (TDHF) frame can be incorporated in the SKKBE by self- 



consistent (i.e. ^-derivable) and conserving approximations 
of the interaction kernel JH). In comparison to (approximate) 
wave function based methods, such as the time-dependent 
Kohn-Sham orbital approach [Toll and its multiconfigura- 
tional variants IfTTlfTZll . NEGF theory does not provide imme- 
diate access to all observables as it is premised on a two-time 
generalization of the reduced (one-particle) density matrix. 
However, all one- and — due to the two-time dependencies — 
also the two-particle quantities can, in principle, be obtained 
at the same level of self-consistency |8] [T3) . A key for ap- 
plying NEGFs to nonequilibrium laser-atom interaction is an 
efficient, well adapted representation of the Green's function. 
In Ref. [ 14 1, we have shown that a suitable concept is given by 
the finite-element discrete variable representation (FE-DVR). 
Here, we extend this work to nonequilibrium. 

It is the goal of this paper to prove feasibility and adequacy 
of two-time NEGF calculations along with spatially extended 
atomic hamiltonians. To this end, we consider N nonrelativis- 
tic electrons which move in a time-dependent potential V(x, t) 
and interact via a pair potential U(x — a/). Sec. |TT] reviews 
the equations of motion for the associated one-particle NEGF 
in FE-DVR representation, where e-e correlations are treated 
in time-dependent second Born approximation (TD2ndB). In 



Sec. Ill we present results for ultraviolet-field (uv) induced 



electron dynamics in case of the helium |[T5l4T9ll and lithium 
hydride model [20 1 including, respectively, two and four elec- 
trons. Thereby, we compare TDHF and TD2ndB to full solu- 
tions of the time-dependent Schrodinger equation (TDSE). 



II. THEORY 

We begin the theory section by briefly describing the grid- 
based representation of the one-particle NEGF G(l, 1') = 

— i(TeV>(l)V ; (!'))> 1 = ( a; i^ <T )' m coordinate space by 
means of a FE-DVR basis — for a more detailed description 
the reader is referred to Refs. [14| and 12H . Thereafter, 
Sec. 
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|ITB 

and Sec. 

memory computation of G(l, 1') using the message passing 



recapitulates the respective equations of motion, 
II C| discusses our concept for parallel distributed- 
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interface (MPI). Atomic units (a.u.) are used throughout, 

TO = |e| = k = 4/T£q = 1, 



A. Representation of G(l, 1') 

We represent the NEGF in one space dimension on an 
interval [0, Xq] which is partitioned into n e finite elements 
[x l , x* +1 ], i G {0, 1, . . . , n e — 1}. The FE boundaries are at the 
same time generalized Gauss-Lobatto (GGL) points x % = x l 



for the one-electron integrals and 



and x 



i+i 



of a GGL quadrature of order n g with ab- 



scissa points x l m and weights w % n 



G {0,1, 



for details see e.g. iFBl and references therein. These points 
and weights define local DVR basis sets |22l . which can be 
superposed such that the total FE-DVR space is spanned by 



TO = 



Xm( x ) 



/A,-i(g) + /o +1 (s) 

K-i + ^o +1 ] 1/2 ' 



(1) 



where f* n {x) are Lobatto shape functions 1221 . All xh( x ) — 
denoted bridge functions — ensure the NEGF being continuous 
in coordinate space, cf. Ref. [14|. Moreover, the basis has 



dimension nj, = n e n g — 1 and is GGL orthonormal |23]. Using 
Eq. Q, the one-particle NEGF G{ 1,1') = 6[t-t') G>(1, 1')+ 
— i) G < (1, 1') in FE-DVR representation is expanded as 

G(xt,x't') = J2T / Xin(x)xi>(x')9% m >(t,t') , (2) 

%%' mm' 

Sw(M')eC, 

where we imply a spin-restricted ansatz, i.e. G does not ex- 
plicitly depend on the spin-projection a. However, the spin- 
degeneracy £ G {1, 2} is incorporated in the self-energies, i.e. 
S(l, 1') = (xt, x't'), the representation of which is as G in 
Eq. (j^ji with time-dependent coefficients ^' mm >(t, t') G C. 

B. Equations of motion 

In FE-DVR representation, the second-quantized form of 
the considered TV-electron Hamiltonian reads as 



■i it, 



(3) 



%v mm 1 

EV^ ?/ *l*2*3*4 [Jl lt|>»*3 if ^2 ^4 

/ ^ li mim2m3m4 I. mi J L mg J m2 7714 ? 

ii...i4 mi ...m 4 



where the operator ([c^]^) creates (annihilates) an electron 
in the local state \ l m (x). In particular, it holds 

= f«' , 4. A"' ,$* (t) (A) 



M mim2 17131714 



= dx dx' xX( x )xlt 3 (x') 



xU(x-x') X % 2 (x)xlr\A X ') 



for the two-electron integrals, where b™7 a — 5u>S mm 
Explicit expressions for the kinetic energy t^ m /, the time- 
dependent potential energy v l m (t) and the interaction matrix 
w" m , are given in Ref. |[T4l . We note, that the equality in the 
last lines of Eqs. Q and Q, respectively, is only true when 
the GGL quadrature rule is being applied. 

Concerning Eq. (|3j, in FE-DVR picture, the equations of mo- 

tionfor 5 « m ,(M') = -i(f c cUt)l4n'(t')V) are the SKKBE 
transformed into the matrix equation 

= S c (t- t')5™ n ' + J (1/ >::',,„„ [g, u] (t, i)9% m > (*. » 

(6) 

where C denotes the complex Schwinger-Keldysh contour and 
additional summation is implied over i and to. The self -energy 
^f mm i[g,u](t,t') depends functionally upon the matrices g 
and u and, in Hartree-Fock and second Born approximation, is 
defined in the Appendix. Due to the simple structure of Eq. (|5]l 
they are easily evaluated in FE-DVR representation. 

The SKKBE (|6|l represents the basic equation to be solved, 
in Sec. Ill for atomic and molecular model systems. Thereby, 
to resolve the correlated dynamics (TD2ndB), it has to be 
propagated into the two-time plane V — [0, tf] x [0, tf] which 
scales quadratically with final propagation time tf. To do so, 
we follow the scheme described in Refs. 



C. Code parallelization and performance 

While the TDHF part of the self-energy, being of first order 
in u, can be absorbed into h" lfh (t) of Eq. (Jy, higher-order 
contributions are connected to e-e correlations and lead to 
a memory kernel [r.h.s. of the SKKBE] which, generally, is 
not smooth. This means, that, during each single time step 
T — » T + St within V, integrals over the complete dynamics 
history have to be performed [25 1. Consequently, for each 
point (t, t') G V an n| -dimensional array has to be allocated 
containing the spatial part of the NEGF. For symmetry reasons, 
the storage of the lesser (greater) correlation function g < (g > ) 
can be constrained to an upper (lower) triangle, cf. the open 
(closed) circles in Fig.[T](a) or (b). Hence, the total amount 
of data evaluates to T> = lQufnl bytes, where n t = tf /St 
and double precision is presumed. In typical one-dimensional 
[atomic or molecular] problems this can easily involve few 
terabytes of data l26l . 

Such memory (RAM) requirements are conceptually beyond 
standard shared memory setups. Therefore, we have developed 
a distributed memory concept [27] based on MPI fl28ll . The 
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FIG. 1: (Color online) Distributed-memory concept in the case of (a) 
two MPI processes, (b) three MPI processes. The different domains 
with dotted, dashed and long-dashed border (thin lines) mark main 
memory for g < (t,t'} allocated by distinct processes (ranks). The 
extent of the memory kernel [the r.h.s. of the SKKBE] for the indicated 
time steps T — > T + 5T (cf. the arrows) is given by the thick dashed 
line, while the required self-energies £<(t, t 1 ) are bordered with a 
thick solid line. Note, that the relevant history is completely known 
to the process that performs the respective time step. In addition, the 
propagation of the greater and lesser Green's function (extending the 
upper and lower triangle) can be performed by a single process. 



main memory allocation for it, t') is embedded as follows: 
Depending on the number of processes available (A4 S j ZC ), the 
two-time plane V is partitioned into a series of different ver- 
tical and horizontal blocks (arrays) attributed to distinct MPI 
processes labeled as -M ra nk in Fig. [TJ More precisely, any 
9mm' (*> * s simultaneously kept in the memory of process 
(tSt^ 1 ) mod A^size + 1 and in the memory of process (t'St -1 ) 
mod Assize + 1. The situation is illustrated in Fig.[TJ(a) for 
the case of two processes, Fig.[TJ(b) refers to Assize = 3. In 
general, we arrive at a (multi-colored) chessboard-like pat- 
tern which has a "primitive cell" of side length (A4 S1ZC — l)St. 
As a consequence, the total NEGF is effectively stored twice, 



i.e. 2? e ff = 21?. However, doing so ascertains that the amount 
of MPI communication is minimized — this is another critical 
issue concerning efficient parallel two-time NEGF calcula- 
tions. In particular, with such memory distribution, no Green's 
function at all needs to be exchanged with other processes ! 
Exemplarily, the area bordered by the thick dashed line in 
Fig. [TJ (a) indicates that, for a selected time step to be per- 
formed by process A'frank = 1 (cf. the arrows), the memory 
[interaction] kernel extends exclusively on local storage do- 
mains. The same holds for diagonal time steps, cf. Fig.[TJ(b). 
Further, note that the propagation of g < and g y can be per- 
formed in sequence but bunched in a single process. As a 
result, only the actual self-energy matrices £5^/ (t, T) and 

E mm' ( T > *') with M' < ^ ( see the areas edged by thick solid 
lines in Fig. [TJ> need to be broadcasted and made known to all 
processes, though they can be computed strictly locally. 

In conclusion, such a distributed memory concept allows 
for an efficient time-stepping algorithm and, at the same time, 
enables large cluster computations. Regarding the scalability, 
NEGF calculations concerning the one-dimensional helium 



atom (cf. Sec. Ill i with up to 512 MPI processes have achieved 
a typical degree of parallelization of more than 95%. Fig. [2] 
shows the corresponding scaling behavior for test calculations 
of different numbers of processes p and length n t . Finally, we 
mention that the memory distribution scheme is easily extended 
to cover also the mixed Green's functions [25 1, which account 
for the self-consistent dynamics of initial correlations. 




2 4 8 16 32 64 128 256 512 
number of processes p 

FIG. 2: (Color online) Speed-up ratio S(p) = Ti/T p (30) for paral- 
lel propagation of G(l, 1') with up to 512 MPI processes for three 
different propagation times tf = ntSt where 8t = 0.025. As physi- 
cal system, we have considered the one-dimensional He model (see 
Sec,|HI|> with nt — 23 FE-DVR basis functions. While the shortest 
calculations (plus symbols, red) have been performed on the "xe" 
nodes of the HLRN batch system [29|, all other calculations have 
been carried out by the "icel" nodes. We note, that the performance 
drop between 4 and 8 processes is caused by internal architecture 
differences. The thin dashed, dotted and dashed-dotted lines indicate 
a degree of parallelization of T = 95% to 99% according to Amdahl's 
law 5r(p). 
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FIG. 3: (Color online) Contour plot (logarithmic) showing the time-dependent one-electron density {h)(x, t) for ID He in Hartree-Fock (solid, 
green) and second Born approximation (dashed, blue). The dashed-dot-dotted (red) lines reveal TDSE solutions, (a) Time evolution of the 
He atom initially prepared in the Hartree-Fock ground state and no uv field being switched on, (b) short-time response of the Hartree-Fock 
initial state to a permanent uv field (£"o = 0.1, ujq = 0.54), (c) permanent uv field as in (b) but with the He atom initially prepared in the 
self-consistent ground state. The contour lines cover a density range from 10 -4 to 0.5 a.u. with four contours in each decimal power between 
10 -4 and 0.1, and contours between 0.1 and 0.5 each 0.05 a.u. 



III. RESULTS AND DISCUSSION 



A. Helium 



To illustrate the method, we consider N electrons moving in 
a one-dimensional model atom or molecule, where the potential 
in Eq. Q plus a permanent laser field in dipole approximation 
is given by 



V(x,t) 



Z m U(x — x Q m ) 

m=l 

E cos(ujQt)9(t)x , 



(7) 



with atomic numbers Z m , M nuclei positions xq . m , a soft- 
Coulomb electron-nucleus or, respectively, electron-electron 
interaction U(x) = [x 2 + l] -1 / 2 , and an electric field strength 
(frequency) Eo (luo). For M = 1, the potential V(x,t) de- 
scribes the core potential of a single atom superposed by the 
laser field, in particular for Z\ = 2 we model the helium atom 
in one spatial dimension (ID He). 



1. Short-time density response 

The neutral ID He model containing two electrons prepared 
in the singlet ground state has total energy —2.2383 (exact). 
In second Born approximation, the ground state energy is 
—2.2334 accounting for 65% of e-e correlations (Hartree-Fock: 
—2.2242), cf. Ref. [ 14|. We now excite this system by uv laser 
light of intensity E Q = 0.1 (3.5 • 10 14 W/cm 2 ) and, from the 
SKKBE resolve the time-dependent one-electron density 
(n)(x, t) = —iG < (xt, x't) within a simulation box of width 
xq = 30 a.u. and up to 49 FE-DVR basis functions. Thereby, 
the frequency luq = 0.54 (84 nm) is chosen such that the 
perturbation is [in all cases] near-resonant to the first excited 
singlet state, cf. Fig. [6] The short-time density response of the 
ID helium atom is shown in Fig. [3] in form of contour plots. 
Fig. [3] (a) shows the dynamics in absence of the laser field, 
E = 0. Here, in TDHF approximation (solid line), the initial 
state remains an eigenstate of H for all times t. However, 
if the Hartree-Fock ground state evolves under influence of 
e-e correlations, see the dashed (dashed-dot-dotted) TD2ndB 
(TDSE) curves, collective oscillations in (n)(x, t) are initiated. 
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FIG. 4: (Color online) Nonlogarithmic plot of the ID He density 
(h)(x, t) for different times t corresponding to Fig.[5](c). The inset 
shows the self-consistent (initial) ground state density at t = 0. Fur- 
ther, the contour lines indicate the absolute of the density difference 
(5h) 1 of TDHF (solid) and TD2ndB (dashed), respectively, to TDSE 
at values of 0.005, 0.01, 0.02, 0.04 and 0.06. 



These small time-dependent features are, in particular, well 
characterized in second Born approximation. 

In Fig.[3](b), laser light is switched on instantaneously at 
t = 0, and, as in Fig. [3] (a), the He atom is prepared in the 
Hartree-Fock ground state. The external field drives the sys- 
tem out of equilibrium leading to strong density deformations 
which reveal that both electrons are oscillating in the uv field. 
Thereby, many time-dependent details in (n)(x, t) [at high den- 
sity around x = 15 and low density] which are not captured 
in TDHF are well resolved in the TD2ndB calculation, see 
e.g. the domains labeled as A, B and C in Fig. [3](b). Also, 
in the TD2ndB and TDSE densities, the oscillatory behavior 
of Fig. [3] (a) is superimposed. Moreover, due to the finite 
simulation box, reflections occur at the interval boundaries. 

Fig. [3] (c) shows the density response after preparation of 
(correlated) self -consistent initial states, i.e. propagation of the 
NEGF in second Born approximation in particular involves 
mixed Green's functions G r (xt, x't 1 ) (G^ (xt, x't')) with t (?) 
on the imaginary branch of the contour C, e.g. (25). These 
can be omitted in case of an uncorrected TDHF calculation. 
Comparing approximate to exact calculations, we see that 
TD2ndB performs fairly well and complies essentially better 
with TDSE than with TDHF. This is observed at high density 
around x = 15 as well as for moderate- and low-density do- 
mains, where (h) (x, t) is spatially more extended in TD2ndB 
approximation. In addition, the oscillatory behavior of (a) and 
(b) vanishes due to the self-consistency. To support the qual- 
ity of TD2ndB, Fig.|4]shows nonlogarithmic snapshots of the 



electronic density close to the nucleus at different times t and 
also quantifies the difference (<5n) 7 = |(n) 7 — (?^)tdse| of 
the approximate to the exact density (solid and dashed contour 
lines). In total, the time-dependent density profile is more 
spread out due to e-e correlations, and, in contrast to TDHF, 
the error (<5«)TD2ndB does not exceed values of 0.05. More- 
over, at specific times during the initial stage of uv excitation, 
at t « 8 and 16 a.u., the deviation to the exact result becomes 
very small. 



2. Time-dependent dipole moment 

>From the NEGF we also can access the time-dependent 
dipole moment (DM) (d)(t) = i J*°dxG < (xt,xt). For the 
helium atom, the DM is displayed in Fig. [5] (a) implying the 
following: (i) To suppress the effect of box boundary reflec- 
tions we have enlarged the simulation box to x = 70 a.u. and 
(ii) have modified the perturbation to a spectrally broad dipole 
kick Eq cos(wf)0(i)x — > E S(t)x with intensity E = 0.01 
(3.5 • 10 12 W/cm 2 ). The latter change allows for a direct 
computation of dipole spectra via Fourier transform of the 
DM time-series. Again, in Fig.|5](a), we compare TDHF and 
TD2ndB to the exact solution of the TDSE. In second Born 
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FIG. 5: (Color online) Time-dependent dipole moment (d)(t) for 
the ID helium atom, (a) TDHF (solid) and TD2ndB approximation 
(dashed-dotted) versus the exact TDSE result (dashed-dot-dotted). 
Note that the time-evolution in case of TD2ndB starts from the 
Hartree-Fock ground state, (b) Comparison of the short-time behavior 
of (d)(t) for different initial states: TD2ndB 2 (dashed-dotted) refers 
to the second Born approximation with Hartree-Fock initial state, 
whereas the initial state for TD2ndB 1 (dashed) has been computed 
self-consistently. 
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FIG. 6: (Color online) Dipole strength (d ) (u>) for ID He. (a) Discrete 
Fourier transform with Hamming window from TDHF (solid) and 
TDSE (dash-dotted) time-series (d)(t) of length tf = 500. The 
first excited state appears at wi = 0.533 (TDSE) and a;^ — 0.549 
(TDHF), the second excited singlet state has frequency UJ2 = 0.671. 
Xi = 0.755 denotes the first ionization threshold. Note, that peaks 
beyond X\ are artificial box states, (b) Auto-regression power spectral 
density (d)s(t) of short TDHF, TD2ndB, and TDSE time-series 
(d)(t) of length tf = 55 computed from Eq. jij with model order 
S = 1000. In TD2ndB approximation, the first excited state appears 

(2) 

at frequency uj{ — 0.537. 



approximation, see the dashed-dotted line, the time-dependent 
DM well approaches the TDSE result with a general shift to a 
larger [main] oscillation period compared to TDHF. Such an 
e-e correlation induced period increase may also be deduced 
by a thorough look at (h)(x, t) during uv excitation of He in 
Fig. [3] (b) and[3](c). Moreover, Fig.[5](b) indicates how the 
evolution of the DM, including e-e correlations, depends on 
the initial state, cf. the dashed and dashed-dotted lines. As ex- 
pected, starting from the self-consistent initial state (TD2ndB 1 ) 
produces a DM which best matches the exact result. In contrast, 
the perturbation of the Hartree-Fock ground state (TD2ndB 2 , 
equivalent to TD2ndB in Fig.[5](a)) produces a different time 
response in the DM amplitude. This, generally, lead [similar to 
the situation in Fig.[3](b)] to additional excitations which are 
not of dipole character. 

A more detailed analysis of the DM is possible in terms 



/ 



of the dipole spectrum (d)(ui) oc f f dt e~ iut {d)(t) with t 
being the final propagation time. However, due to the two- 
time structure of the SKKBE ([6]l, typical TD2ndB calcula- 
tions are limited to relatively short times of approximately 
< 100 a.u. even if the parallel algorithm of Sec. Ill C : 



t 



is be- 



ing applied. Consequently, a discrete Fourier transform (DFT) 
reveals poor resolution, cf. the dashed line in Fig. [6] (a). In 
the case of TDHF (solid line) and TDSE (dashed-dotted line) 
much better resolution can be achieved for (d ) (w) such that, 
in Fig. [6] (a), one recovers the series of one-electron excitations 



Wi (excited singlet states of ID He) which converges towards 
the first ionization threshold X\ = 0.755. The peaks beyond T\ 
indicate transitions to eigenstates of the finite simulation box. 
Overall, there is a shift concerning the first excited singlet state 
while energetically higher excitations seem to be less affected 
by e-e correlations. Concentrating on the first excited state of 
frequency ui\, we observe that it is overestimated by TDHF 
(u>i — 0.549) — the exact value is 0.533. Thus, it is inter- 

(2) 

esting to analyze the behavior of TD2ndB (oj\ '). However, 
as usual Fourier transform of the DM time series is not well 
applicable, we have to rely on other methods which are not 
limited by the total length of the signal, e.g. harmonic inver- 
sion (HI) l32ll . Nevertheless, the best results have not been 
obtained by HI techniques, but by using an auto-regression 
(AR) spectral analysis [33| where an auto-correlation model 
of order S is applied to the time-dependent dipole moment: 



(d)(t) = (d ) + J2 c s{d)(t - a At) , At = 



f 



(8) 



with time average (do) = dt (d ) (t) and AR coefficients c s 
which, e.g., follow from the Yule-Walker equations [33 1. The 
AR coefficients then give direct access to the power spectral 
density (PSD) via 



(d)s(u) oc 



s 

1 — c s exp(— isw<5i/2) 



(9) 



which is periodic in lo$ = 47rAi~ 1 = AirS/tf. Fig. |6](b) 
summarizes the AR PSD result for DM time-series of fixed 
length tf = 55 OTI . Concerning the first excited state, TDHF 
(solid line) and TDSE (dashed-dot-dotted line) lead to well pro- 
nounced peaks which exactly equal the DFT result of Fig.[6](a). 
In the same manner, we are now able to explore the TD2ndB 



data and find lo^ = 0.537, cf. the dashed line in Fig. [6] (b). 



This frequency corresponds to a remarkable improvement of 
75% of the deviation Hartree-Fock to exact. In addition, in 
(d)s(w), also the second excited singlet state is resolved. The 
general shift towards smaller frequencies 0J2 is, thereby, caused 
by the shortness of the time series. In summary, the AR spec- 
tral analysis, being superior to DFT, allows for an appropriate 
characterization of short (d)(t) time series. In particular, we 
conclude that, concerning dynamic properties of the ID He 
model, TD2ndB shows correlation induced features to a sub- 
stantial level, though, in the present analysis, it was not possible 
to show how well two-electron excitations fi"2ll — emerging be- 
tween 2i and I2 — 2.2383 and being absent in TDHF — are 
resolved in TD2ndB approximation. This lack is only due to 
the limited propagation time, and, indeed, we expect these 
features of e-e correlations to be included on a similar level. 



B. Lithium hydride 

Using two-time NEGFs, it is possible to describe even more 
complex systems than He in nonequilibrium. As a proof of 
principle, we consider the response of the neutral molecule 




5 10 15 20 25 30 

position x (a.u.) 



FIG. 7: (Color online) Time-dependent one-electron density {h)(x, t) 
(plotted nonlogarithmically) for the response of LiH to an xuv field 
with intensity Eo = 0.75 and frequency cjo = 1.3. While panel (a) 
indicates the [exact] TDSE result, panel (b) and (c), respectively, refer 
to self-consistent TD2ndB and TDHF calculations. Arrows in (c) 
indicate particularly strong deviations of TDHF compared to TDSE. 



lithium hydride (LiH) to an xuv field of intensity Eq = 0.75 
(2.0 • 10 16 W/cm 2 ) and frequency uj = 1.3 (photon energy 35.5 
eV). LiH is a heteronuclear molecule (M = 2, Z\ = 3, and 
Z 2 = 1 in Eq. (|7]i) including four electrons and complies with 
the restricted NEGF ansatz Q when prepared as singlet state, 
i.e. with closed-shell spin configuration. Prior to excitation, we 
have fixed the Li-H bond lengths — xo,2 — ^o.i t m terms 
of the Born-Oppenheimer separation of electronic and nuclear 
motion] to the self-consistent values as obtained by scanning 
the potential energy surface lfT4ll : db — 3.386 (Hartree-Fock), 
dh = 3.505 (second Born), db ~ 3.6 (exact). The resulting 
ground state densities (n) (x, 0) are indicated in Fig. [8] cf. the 
lines at the bottom of the plot (t = 0). Computing the LiH 
time evolution from the TDSE in four spatial coordinates is, in 
contrast to helium, already ambitious but nevertheless feasible. 
In comparison, going from ID He to the LiH model solving the 
SKKBE for the one-particle NEGF is conceptually and compu- 
tationally very simple as, with the increase of particle number, 
only the normalization of the Green's function is concerned. 
Fig.[7]reveals the picture of the xuv induced electron dynamics 
over more than three laser cycles as obtained from TDSE (top), 
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FIG. 8: (Color online) One-electron density (h)(x,t) at different 
times t for the response of LiH to an xuv field (with parameters as in 
Fig.|7}. The TDHF (solid) and TD2ndB 1 result (dashed) include self- 
consistent initial states (with different but fixed db) while TD2ndB 2 
(dotted) starts from the Hartree-Fock ground state [only shown for 
t < 10)]. Note, that the self-consistent solution of the TDSE (dashed- 
dot-dotted) is not fully accurate due to lack of spatial resolution — 
compare with the exact ground state (dashed-dotted line) for t = 0. 
Moreover, the gray curves show the time-dependent expectation value 
(x) (t) of the electron position. 



time-dependent [self-consistent] second Born approximation 
(center) and time-dependent Hartree-Fock approximation (bot- 
tom). Thereby, we have used an FE-DVR basis of size rib = 71 
and a simulation box of xq = 40. As general trend, we observe, 
in (n}(x,t), that TD2ndB is considerably superior to TDHF. 
In particular, the time-dependent density between the Li and 
the H atom behaves completely different than TDSE in the 
case of TDHF, see e.g. the density around x — 20 for t = 9 
and t = 13 indicated by arrows. This is essentially corrected 
in TD2ndB approximation. 

For specific times t, Fig. [8] allows for a direct and more 
detailed comparison of the exact and approximate one-electron 
densities close to the atomic nuclei. As in Fig.pl TD2ndB 1 
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(dashed-dotted line) means that the initial state is uncorrelated 
[thus on the Hartree-Fock level], while TD2ndB 2 (dashed line) 
refers to self-consistent ground state preparation. Following the 
time-evolution, TD2ndB 2 leads to the correct changes towards 
the TDSE solution including strong depopulation of electrons 
around the hydrogen atom for t > 10. Also, the TD2ndB 2 den- 
sity is more smooth whereas TDHF shows additional maxima 
being absent in, both, TD2ndB approximation and the exact 
result (dashed-dot-dotted line), see e.g. the density for t = 7.5 
at x 20. In addition, we observe that starting from the corre- 
lated (TD2ndB 2 ) or uncorrelated ground state (TD2ndB 1 ) does 
not make much density difference for later propagation times 
t > 5. Averaging over coordinate space allows us, moreover, 
to trace the time-dependent expectation value of the electron 
position (x) (t) which, being initially close to the Li nucleus, 
starts to oscillate with the xuv field, cf. the gray lines in Fig. [8] 
Thereby, TD2ndB, again, gives consistent results, although, 
for later times, deviations to TDSE increase. To this end, we 
note that the solution of the full four-particle TDSE could have 
been performed only with restricted spatial resolution whereas 
the NEGF results are converged. Thus, comparisons should 
be drawn carefully. Nevertheless, for confirming the main 
trends [when e-e correlations are being included in the NEGF 
calculation], the TDSE result is fully adequate. 



atoms or molecules [where exact calculations are impossible] 
affirms the power of the present NEGF approach. 

As the TD2ndB results show important time-dependent fea- 
tures, we believe, that the present analysis will stimulate and 
support further investigations of strong-field laser-atom inter- 
actions based on the two-time NEGF approach. Moreover, due 



to the distributed memory concept (Sec. II C i which is capable 
to organize the huge demand of RAM involved in the memory 
kernel, we are now able to perform the two-time propagation 
of the SKKBEs in an optimized and efficient way. This issue, 
in future work, will enable more extended calculations with in- 
creased FE-DVR basis size {rib) and longer propagation times 
(t f) and, hence, fortifies a quantum kinetic approach to atoms 
and molecules and their interaction with (strong) laser fields. 
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IV. CONCLUSION 



Appendix: Self-energies 



The (x)uv laser-induced, short-time electron dynamics in the 
ID helium and lithium hydride model have been investigated 
using a NEGF approach. To this end an efficient grid-based 
approach of Ref. [ 14] was extended to nonequilibrium situa- 
tions. Starting from initial states with different degrees of e-e 
correlation and self-consistency, it has been shown that NEGFs 
are well applicable to resolve the time dependence and that 
dynamic correlation effects are approximated to a substantial 
and valuable level already in the case of TD2ndB. This is ob- 
tained from (i) analyzing the few-cycle (one-electron) density 
response of He and LiH to strong laser fields and (ii) from 
exploring spectral properties such as the singly-excited spec- 
trum of ID He, cf. Sec. [Hi] Overall, results at TD2ndB level 
are found to be closer to exact solutions of the TDSE than to 
TDHF. Thereby, the present restriction of NEGF calculations 
to examples of two- and four-electron systems stems from the 
availability of exact reference data which allows for a proof 
of principle. Exact solutions of the TDSE for N > 4 are es- 
sentially at the limit concerning spatial resolution whereas a 
solution of the SKKBEs is not explicitly sensitive to the parti- 
cle number. This potential of direct extension to larger model 



For a complete discussion of the matrix equation (|6]l in 
Sec. II B we, here, give the definitions of the self-energies: (i) 
In time-dependent Hartree-Fock approximation (TDHF), the 
self-energy is given by T,^[g, u] = 8c(t — t')T^ F [g, u] with 

— -\fS ili2 v ilia a* 3 ' 3 (t /+) 

+ iu l2H a i2il it t + ) 

' 1 u, m 2 "il i)m 2 m\ \ L i L J ' 

(ii) in time-dependent second Born approximation (TD2ndB), 

it is E € [s,u] = 5 C (t - i')£f F [g, u] + Zf ldB {g,u}, where 

eK"' [»,«](*,*') 
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